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Abstract 

Inclusive single-particle production cross sections have been calculated in- 
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are compared with experimental data from the CERN SppS Collider and the 
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1 Introduction 



The inclusive production of single hadrons in hadronic collisions has long been recognized 
as an important process for testing the QCD-improved parton model (for reviews on the 
subject see, for instance, ref. [IH). The leading-order (LO) QCD formalism consists of 
using the tree-level results for the hard cross sections, the one-loop expression for the 
running coupling constant, and parton densities (structure functions) and fragmentation 
functions generated by one-loop evolution kernels. This theory provides a rather consistent 
description of many large- momentum-transfer processes 0. The tree- level parton-parton 
scattering cross sections that enter this formalism have been calculated many years ago 
and an extensive phenomenology has been developed and confronted with experimental 
data 

First results for next-to-leading order (NLO) corrections to hard-scattering processes 
have been obtained by Ellis et al. in 1980. However, these authors considered only 
subprocesses with quarks in the initial state. Using quark distribution functions defined 
in deep inelastic scattering and fragmentation functions defined in one-particle inclusive 
e"'"e~ annihilation, they found large corrections. Meanwhile, a complete NLO calculation 
involving all parton-parton processes contributing to hadron-hadron collisions has been 
presented by Aversa at al. [Q. 

These results are requisite in order to perform a consistent calculation of single-hadron 
production in the NLO formalism. This involves NLO hard-scattering cross sections with 
two-loop as, two-loop-evolved structure functions and fragmentation functions. Only 
under these circumstances are all four elements unambiguously defined. Using NLO results 
for hard-scattering cross sections, structure and fragmentation functions in combination, 
should reduce substantially the scale dependence of the predictions for inclusive one- 
particle cross sections, since the variation in one of these factors is compensated by the 
others up to still missing higher-order terms. To guarantee this compensation, a consistent 
choice of renormalization and factorization schemes and scales is necessary in all four 
elements entering the analysis. As mentioned above, NLO results for all parton-parton 
subprocesses are available. 

As is well known, structure and fragmentation functions cannot be calculated yet 
from the QCD Lagrangian. For their construction one has to resort to experimental data, 
in particular from total inclusive deep-inelastic lepton-nucleon scattering in the case of 
the structure functions and from one-particle inclusive e^e~ annihilation in the case of 
the fragmentation functions. As for structure functions, such analyses based on NLO 
evolution and all existing data have been carried out by several groups, even including 
processes other than deep- inelastic lepton-nucleon scattering |0, 0; for a recent review see 

i- 

Unfortunately, such analyses do not exist for the fragmentation functions of quarks 
and gluons. Although the relevant NLO theory is well known |TD[ and inclusive single- 



particle distributions of various hadrons have been measured in e^e~ annihilation as well 
as in semi-inclusive deep-inelastic lepton-nucleon scattering, such a NLO analysis has not 
been tackled yet. All the existing information on fragmentation functions is based on 
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LO evolution equations and also on rather old data 0, |TT|. Although this represents a 
serious limitation to any NLO prediction of one-particle production in hadron-hadron, 
photon-hadron, or lepton-hadron reactions, we consider it appropriate to make the first 
step and to combine the available results for the NLO parton-parton hard-scattering cross 
sections and the NLO parton distribution functions to obtain more reliable predictions as 
compared with the LO theory considered in the past. 

Inclusive single-hadron production provides both advantages and disadvantages as 
compared to jet production. A distinct hadron constitutes a final state which is better 
defined than a jet event, since the latter delicately depends on the details of the jet- 
finding algorithm. Moreover, the definition of a jet becomes more problematic at low 
energy. On the other hand, the fragmentation of the produced parton renders the px 
distributions quite steep and reduces therefore the range experimentally accessible in the 
case of single-hadron production. Accordingly, the experimental errors are large in the 
high-pT range, where the cross sections become considerably small, while the bulk of the 
data is accumulated at lower values of p-p. 

In this paper we shall present cross sections for single-charged-hadron and single-7r° 
production in pp reactions. Similar predictions for 7p and low-Q^ ep processes will be 
considered in a separate communication. Experimental data for charged one-particle and 
tt" production at relatively low centre-of-mass (CM) energies come from the CERN ISR 



pp Collider [12|. At higher energies, data from the UAl [13] and UA2 [14| Collaborations 
at the SppS Collider and from the CDF Collaboration [|T^ at the Tevatron are available. 
Since the QCD-improved parton model should work best for the pp collider energy range, 
we shall compare our results with data for the higher energies only. Specifically, we choose 
data at ^/s = 540 and 630 GeV from UA2 and CDF, respectively and at 1.8 TeV from 
CDF. 

The outline of this paper is as follows. In Sect. 2, we give a short introduction of 
the formalism and fix our notation and input. Results of our calculation are presented in 
Sects. 3 and 4, where we also compare with data from pp colliders. Sect. 5 is reserved for 
a discussion of the results and some concluding remarks. 
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2 Formalism, Notation, and Input 



In the LO QCD-improved parton model the single-particle cross section can be ex- 
pressed as a convolution of the LO parton-parton scattering cross sections with the scale- 
dependent structure and fragmentation functions. We fix the momenta of the hadron 
reaction by 

/il(pi) +/i2(p2) /i3(P3) +X, (1) 

where hi and h2 are the incoming hadrons and is the single hadron in the final state. 
In LO the cross section is given by 

Y.ldxJ dx, [ ^ F^^{x,^M')F^\x,,M')D^^ (x,,Mf) pl^-^^^ (2) 

where the partonic momenta ki, kj and ki are expressed in terms of pi, p2 and ps as: 

ki = xipi, kj = X2P2, h = P3/X3 . 

The indices i,j, I run over gluons and Np flavours of quarks. We assume Np = 4 through- 
out our calculation and neglect the influence of the charm-quark threshold. F^{x,M'^) 
and Dj^{x,Mj) are the usual structure and fragmentation functions appropriate to par- 
tons of type i or / inside hadron h, which depend on the scales M and Mj to be specified 
later. Finally, d^a^^.f^.^i^^ characterizes the process i + j^l + X in O (a;g(/i^)). 

Thus, we have three scales altogether: M, the factorization scale of the parton dis- 
tributions F^, Mf, the factorization scale of the fragmentation functions D^, and fi, the 
renormalization scale of the QCD coupling constant as. In LO all three scales are left 
undefined. Usually they are linked to the transverse momentum of the produced 
hadron h^. Since this choice of scale is rather arbitrary, there is a large uncertainty of the 
theoretical predictions at the LO level. 



In NLO the inclusive cross section is written a^ 
^3 



d^p, 



'3 



Y: I dx.J dx2j '^Fl'^{xi,M^)F^\x2.M^)D^- {x„ Mf 



X 



-3 

- I \5{l-w) + ^— Kk^k.^k, [s, v,w;iJ. ,M ,Mf) 



(3) 



where v and w can be expressed in term of the partonic variables s, t and u as: 

v = l + t/s, w = —u/{s + t). 
tin the following, we largely follow the notation of Aversa et al. M. 
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In turn, s, t and u are related to the hadronic variables S 
and U = {p2 — VzY by: 



X1X2S . 



t 



Xi 



T . 



u 



X2 



U 



The factors K, 



M- 



are the NLO corrections. 



The calculation of the correction terms K^.^.^k^ is fully described in 0. It is based 
on the results of Ellis and Sexton for the matrix elements squared for all (2 2) 
and (2 3) parton scattering subprocesses calculated to O (ag). The Kk^kj^ki functions 
are plagued by infrared and collinear divergences, which appear as poles in the parameter 
e = (4 — n)/2, with n being the dimensionality of space-time. The cancellation of the 
infrared singularities, which are oc is explained in [0. The left-over 1/e singularities 
associated with ingoing parton lines are absorbed into the bare distribution functions 
Fj^(x), which renders them dependent on the scale M. Similarly, the 1/e singularities 
from outgoing parton lines are absorbed into the bare fragmentation functions D^[x), so 
that they become dependent on the scale Mf. 



After all these cancellations and absorptions have been done, the single-hadron in- 

ki being given by 



elusive cross section has the form (3) with the correction terms Kkikj 



Kk^k.^ki {s,v,w;ij\M\Mf 

1 / dalk,^ki{s,v) 
V \ dv 




s ~ s s \ 

ci + ^iln^ + 5iln^ + ciln— j 6{l-w) 



+ (^2 + 52ln^ + 52ln^ 



1 



+ C3 



1 — If / + 



ln(l — w) 

1 — w 



(4) 



where the Ki 



[s,v,w) functions are regular in the limit w 1. The coefficients q, 
Cj, Cj, and q may be found in for the various subprocesses. The K'^^^,^^^^ functions are 
given by rather lengthy expressions and have been communicated to us by the authors of 



The form of the various coefficients and the K'j^.f^.^i^^ functions in (^) depend on the 
choice of the finite corrections fij{x) and dij{x) = q, g) to structure and fragmentation 
functions. These finite corrections have to be adjusted according to the selected factoriza- 
tion scheme. The authors of advocate a scheme where all fij and dij are kept different 
from zero in the form specified in [Q. However, this alternative reaches beyond the scope 
of present structure-function analyses, since it affects higher orders relative to the NLO 
level, and we shall not follow it. The popular DIS and MS schemes require fij = dij = 0, 
except for fqq and dqq, which are finite and different for the two schemes. The expressions 
for fqq and dqq listed in correspond to structure and fragmentation functions in the MS 
scheme. Since fragmentation functions have been determined only to LO, the MS scheme 
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is appropriate here {dqq 7^ 0). As for structure functions, results are available in both MS 
and DIS schemes. Depending on the choice of scheme, the appropriate forms of fgg must 
be selected in the NLO corrections to the parton-level cross sections. 

As already mentioned in the introduction, we shall employ the fragmentation functions 



by Baier et al. B and Anselmino et al. which appears to be the best choice for the 



time being. Reference []TT| provides us with functions for the fragmentation of u and d 
quarks into vr^ and as a function of the scale Mf. They are obtained from fits to 
EMC deep-inelastic lepton scattering data. We adopt the fragmentation functions of the 
gluon into tt^'*^ from |Q, with D^^ = = . These, as well as the u and d quark 
fragmentation functions have a reference scale Qq = 5 GeV. Following [|l7l, we identify the 
fragmentation function of the s quark into a given pion with the one of the first-family 
quark that is not a ground-state constituent of the respective pion, e.g., D'^^ = -DJ^, etc. 
We proceed similarly in the case of c quarks. The fragmentation functions into kaons 



emerge from those into pions through multiplication by heuristic reduction factors [11 



The proton structure function parametrizations used in our calculcations are taken 



from the package by K. Charchula []TB|. For the central results and, in particular, for the 
comparison with experimental data we use the Morfin-Tung set MT-Bl 0, either in the 
DIS or in the MS version. 

In the next sections, we shall present our results in two steps. In the first step, we 
assume D'l{xs, Mj) = 6{1 — x^) for all hadrons h and partons i. In this way, we study 
the cross sections for the production of single partons. Although these cross sections 
cannot be compared with experimental data — a fragmentation function ~ 6{1 — X3) is 
a reasonable approximation only for heavy quarks — , they are simpler to calculate and 
reveal several features and dependences that are also relevant for the complete single- 
hadron cross sections considered in the second step. We emphasize that these cross 
sections must not be compared with experimental jet production cross sections, since 
they lack a recombination algorithm indispensable for defining jets. 
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3 Results: Single-Parton Inclusive Cross Section 



In this section we approximate all fragmentation functions by 6 functions. This leads to 
single-parton inclusive cross sections. For a first orientation, we disentangle the various 
parton-parton channels that contribute in LO to the hard-scattering cross section, viz. 

1) 99 99, 

2) 99 QQ, 

3) 9(1 9(1, (19 and qg qg, gq, 

4) qq^99, 

5) qq qq, 

6) qq qq, qq, 

7) qq q'q', 

8) qq' qq', q'q, 

where q and q' denote different quark flavours. In Fig. ||, this is done for da/dp^ at 
y/s = 630 GeV using the LO set SL of the Morfin-Tung structure functions 0. Here and 
in the following figure we choose /i = M = px- We observe that at small Pt the two 
channels 1) and 3) dominate by more than one order of magnitude, whereas at larger 
(Pt > 30 GeV) also the pure quark channels 6) and 8) contribute significantly. Similar 
observations have been made in the past It is clear that the magnitude of the cross 
sections with gluonic initial states are highly sensitive to the gluon content of the proton 
(antiproton), whose small- a; behaviour is still not well known. Of course, the relative 
size of the various channels in Fig. |1] is subject to change when realistic fragmentation 
functions are taken into account. 

Next, we study how our results are infiuenced when different sets of structure functions 
are used. Toward this end, we calculate d^a/dyd'^pT in the Born approximation at ^/s = 
630 GeV and y = as a function of using in turn the structure-function sets MT-Bl 
(DIS), MT-Bl (MS) and MRSn-DO (MS), MRSn-D- (MS) which are all NLO 
structure functions. For the QCD coupling as we employ the usual one-loop formula. In 
Fig. 1^, we show these cross sections normalized with respect to the value obtained when 
the MT-SL set is used. Strictly speaking, only the latter should be combined with the 
LO hard- scattering cross sections because it is a genuine LO fit; all the others correspond 
to NLO analyses of present data. 

We observe that the variation caused by replacing the LO set by the various NLO sets 
does not exceed some 30% for px > 5 GeV. The DIS and MS versions of MT-Bl differ by 
only 10% relative to each other. We may consider this as an upper bound on the scheme 
dependence, since the partial compensation that is expected to take place at NLO due 
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Figure 1: pt distributions of the various LO channels at -y/s = 630 GeV evaluated with the LO 
set SL of the Morfin-Tung structure functions. 
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Figure 2: Influence of typical NLO structure functions on the LO calculation of d^a/dy (Ppx 
at y/s = 630 GeV and y = 0- The curves are normalized to the LO calculation with the MT-SL 
set. 



to scheme-dependent terms is not in effect here. The MRSn-DO and MRSn-D— fit also 
recent NMC and CCFR data which extend to smaller values of x than previously 



tested. The results differ significantly at low pt- In particular, the DO and D— sets differ 
drastically from each other. This may be understood by recalling that the gluon-initiated 
subprocesses dominate at low p^', see Fig. [l|. The variation is considerably reduced at 
larger p^ (pt > 30 GeV), where also the channels (6) and (8) contribute significantly. 
These features are expected to survive at NLO. 

To investigate the influence of the higher-order terms in the hard-scattering cross 
sections, we calculate d^a/dy d'^p^ at ^/s = 630 GeV and ?/ = as a function of 2pt/v^ 
with and without these terms and plot the ratio in Fig. ^. In both numerator and 
denominator we use the same structure functions, namely the DIS or MS versions of MT- 
Bl, and the two-loop expression for ag. The renormalization scale A of as is adjusted 
according to the structure functions. We do this for the three choices fi = M = Mj = 
Pt/2, Pt, and 2pT- From Fig. ^ we see that this ratio is independent of p^, except for 
small Pt- Its value ranges between 1.3 and 2.1, which demonstrates the importance of the 
NLO corrections to the parton-level cross sections. The scheme dependence DIS versus 
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Figure 3: d?a/dy cPpx in NLO normalized to the corresponding calculation with the NLO terms 
in the hard-scattering cross sections omitted, for ^/s = 630 GeV, y = 0, and scales set to pt/2, 
Pt, and 2pT- DIS and MS results are compared. 

MS is negligible. Note that, in contrast to Fig. 2, the same structure functions are used 
in both numerator and denominator, so that variations due to changes of scheme or set 
largely cancel out. On the other hand, we observe a strong scale dependence, which is 
mainly due to the denominator as will become apparent below. 

The scale dependence of numerator and denominator of the ratio considered in Fig. 3 
is explicitly studied in Fig. ^ and |b for px fixed at 31.5 GeV and 94.5 GeV, respectively. 
We set the three scales /i, M and Mf equal to and we let ^ vary between 0.125 and 
8. The higher values of ^ in Fig. allow for scales exceeding the CM energy. For 
the results obtained in these figures, the DIS scheme is adopted in connection with the 
appropriate MT-Bl set. Obviously, the LO cross section (dashed line) exhibits a stronger 
dependence on the scale factor ^ as compared to the NLO result (solid line). The more 
rapid decrease of the LO cross section for increasing ^ explains the mild increase of the 
ratio plotted in Fig. ^. The NLO cross section is particularly insensitive to ^ for ^ > 1. 
This holds true independently of px, although the absolute cross section strongly varies 
with Pt] the cross sections shown in Figs. ^ and b differ by four orders of magnitude. 
The scale variation of the NLO cross section is well bounded in contrast to the LO cross 
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Figure 4: Scale dependence of d^a /dy d?pT for y/s = 630 GeV, y = 0, and a) px = 31.5 GeV, 
b) Vt = 94.5 GeV. The solid and dot-dashed lines represent the NLO results with variable and 
fixed Mf, respectively; the dashed lines represent the LO result with NLO structure functions 
and as. 
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section, which decreases monotonically with increasing ^. However, one should keep in 
mind that the Mf dependence of the NLO result is not compensated as long as (5-function- 
type fragmentation is assumed. In this somewhat pathological situation it would be more 
appropriate not to vary Mf. For comparison, we have included the respective results in 
Fig. ^ also (dot-dashed curves). 



11 



4 Results: Single-Hadron Inclusive Cross Section. 



In this section we shall present our results for the inclusive production of charged hadrons 
and of single vr^'s. For the calculation of these cross sections we have employed the frag- 
mentation functions already mentioned in Sect. 2. For the comparison with experimental 
data we have chosen the structure function MT-Bl in the DIS scheme. Collider data for 
charged hadrons {h = {h^ + h')/2) come from the UA2 and CDF Collaborations |15 



In Fig. 1^ we compare the charged-hadron inclusive cross section as measured by UA2 at 
^/s = 540 GeV, for 1.0 < y < 1.8, with the results we obtain when the three different scales 
/i, M and Mf are simultaneously set equal to Pt/2, Pt and 2pT- This scale dependence is 
rather moderate and we shall discuss later on how it compares with the scale dependence 
in the LO cross section. Up to this uncertainty, this is an absolute prediction of the cross 
section. We observe that the prediction with scales equal to px is in adequate overall 
agreement with the data for up to 10 GeV. The agreement is less satisfactory at small 
Pt, where higher order corrections may be appreciable. In addition, somewhere in this 
region non-perturbative effects come in. 

Similar conclusions can be drawn for the comparisons shown in Figs. ^ and |^ with 
CDF data at ^/s = 630 GeV and = 1.8 TeV, respectively. In both cases, the rapidity 
is averaged over the interval —l<y<l. The 1.8 TeV data at large pt are again predicted 
best with scales equal to px- The same seems to be true for the 630 GeV data, although 
the Pt range is in this case too small to allow a firm statement. 

Looking at the energy dependence, we note that y/s increases by a factor 3 as we go 
from Figs. |^,^ to Fig. 0. We observe that at small px, below 3 GeV, the cross section is in- 
dependent of ^/s, whereas in the large-px range the cross section increases with increasing 
energy, i.e. it shows the characteristic flattening of the px distribution originating from 
the hard-scattering cross section. This behaviour of the data is borne out by the theory 
quite well although at small pt we would not expect such a good agreement. 

In order to get an idea of the size of the NLO correction in the hard scattering terms, 
we compare in Figs. ^, ^ and |TD| the Born cross section with the NLO one for y/s = 



540, 630 GeV and 1.8 TeV. Here, again, the cross sections are averaged over the rapidity 
ranges in the same way as for the respective experimental data. In the Born cross section 
the structure functions as well as the value of as are the same as in the NLO cross section. 
The effect of the NLO terms leads to an average K factor of approximately 1.5 for the 
choice of scales /i = Mf = M = p^- For comparison, we show also the experimental data. 

It has to be emphasized that this K factor is only a measure of the NLO corrections to 
the hard scattering. Usually one takes the attitude of comparing the NLO cross section 
to the LO one where the hard scattering, the structure functions and the strong coupling 
constant as are evaluated in lowest order. Such a K factor, denoted as full/ "born", is 
presented in Fig. p!T| a for the three scales pt/2,pt and 2pT, CM energy i/i = 1.8 TeV 
and \y\ < 1.0. The parametrization MT-SL is used for the structure functions entering in 
the LO cross section. The value of 144 MeV for the QCD parameter A required by this 
parametrization is extended also to the as coupling constant and to the fragmentation 
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pp -» h + anything 




3 6 9 

Pt (GeV) 

Figure 5: Inclusive charged-hadron production cross section for ^/s = 540 GeV, averaged over 
the rapidity range 1.0 < y < 1.8. The short-dashed, solid and long-dashed lines correspond to 
the full NLO predictions for scales fi, M and Mf simultaneously set equal to Pt/2, and 2pT- 
For comparison, the UA2 data taken at the same CM energy and in the same rapidity range are 
also shown. 
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Figure 6: Same as in Fig. |5| for ^/s = 630 GeV and rapidity range \y\ < 1.0. The theoretical 
resuhs are here compared with data obtained by the CDF Collaboration. 
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Figure 7: Same as in Fig. | for ^/s = 1.8 TeV. 
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pp ^ h + anything 
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Figure 8: The experimental data shown in Fig. ^ are here compared to the NLO theoretical 
result (upper curve) and the LO one (lower curve) as obtained for structure functions and as 
kept at the NLO level. 
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Figure 9: Same as in Fig. ^ for the data shown in Fig 
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Figure 10: Same as in Fig. ^ for the data shown in Fig. 0. 
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Figure 11: Ratios of the NLO charged-hadron production cross section over a) the LO one and 
b) the LO one evaluated with NLO structure functions and qs, for ^/s = 1.8 TeV, rapidity range 
\y\ < 1.0 and scales equal to pt/2, pt and 2pT- 



functions. We observe in Fig. [TT| a that the K factor obtained in this way is ~ 1 and almost 
PT-independent for > 5 GeV if the three scales fi, M and Mj are set equal to p^. For 
the other two scales, Pt/2 and 2pT, the K factor has a somewhat stronger dependence 
and reaches 0.9 and 1.4, respectively, for the higher values of p^, where experimental data 
are available. 

The K factor with structure functions and as kept at the NLO level, referred to above 



{full/horn), is shown in Fig. |TT|b also for the same choices of scales. It is larger than the 
one shown in Fig. |Tl|a, ranging between 1.2 and 1.8 at the upper end of the px-range 
studied here, and exhibits also a stronger scale and dependence. We have verified that 
similar features are present also at y/s = 540 and 630 GeV. K factors between 1 and 2 are 
obtained in these cases for scales varying between pt/2 and 2pT, for 5 GeV< < 10 GeV. 
A weaker px dependence is observed at these lower CM energies: the ratios full/born are 
in these cases essentially fiat for px > 5 GeV as also appears in Fig. ^, where (5-function 
fragmentation functions are assumed. 

It is clear from these figures that if one goes from a bona fide LO theory to a NLO 
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combination of hard-scattering cross sections, structure functions and as, a strong com- 
pensation between these NLO effects takes place. On the other hand, less compensation is 
present when the NLO effects are contained only in the hard scattering. The importance 
of including the NLO corrections consistently in all components, as already emphasized 
by many authors (see for example f^), is therefore quite evident. 

To obtain more information on the scale dependence of the inclusive cross section, we 
have calculated it as a function of the scale factor ^ defined above for = 5, 10 and 
30 GeV, for CM energy ^/s = 630 GeV and rapidity range \y\ < 1.0. Results are shown in 



Fig. [T^ for the LO cross section (with LO structure function and as), for the NLO cross 
section and for a LO cross section with the same structure functions and as used in the 
NLO calculation. 

We see quite clearly that the NLO cross section has much less scale dependence than 
the two LO ones. The NLO cross section develops a plateau at small values of ^, for not- 
too-small pt, i.e. pt > 5 GeV. The large turn-over of this plateau, however, is presumably 
an artifact due to the use of LO fragmentation functions. For all p^, both LO curves show 
a monotonic dependence on ^. Between ^ = 1 and ^ = 8 the decrease of these two curves 
is steeper than the one shown by the NLO cross section. We also observe that the LO 
cross section (dotted curve) equals the NLO cross section for all three px at ^ ^ 0.7, 
whereas at larger values of ^ the two curves are quite well separated. This scale, where 
the true LO and the NLO results match, is advocated as the best possible choice of scales 



by the method of Fastest Apparent Convergence [22 



At this point, we should compare Fig. |T2| for p^ = 30 GeV with Fig. where the 
scale dependence has been investigated at y/s = 630 GeV and pr = 31.5 GeV in the case 
of 5-function fragmentation. We notice that, apart from the absolute size, the LO results 
exhibit roughly the same qualitative features, which can be attributed to a rather weak 
scale dependence of the fragmentation functions. However, the situation changes appre- 
ciably when NLO corrections are included. Then, in the case of realistic fragmentation 
the cross section falls off less steeply at large ^, due to compensations between the Mf 
dependence of the fragmentation functions and the Mf dependence of the hard-scattering 
cross sections. Finally, the drop-off at small ^, when present, is also more pronounced in 
Fig. 0. 

So far we have concentrated on the discussion of pt distributions with rapidity inte- 
grated over a finite region. In Fig. |1^ we show the differential cross section d^a/cPpTdy 



as a function of the rapidity y for CM energy ^/s = 630 GeV and three different values of 
Pt'- Pt = 5, 10 and 30 GeV. We observe that the NLO corrections to the hard-scattering 
cross section do not change the shape of this distribution. The full NLO cross section and 
the LO one both exhibit the same symmetric behaviour in y. 

At last, we show the results we obtain on inclusive tt^ production at ^/s = 540 GeV 
for 1.0 < y < 1.8 (Fig. 14) and \y\ < 0.6, compared with the UA2 data taken at the 



same CM energy and in the same regions of rapidity (Fig. ^). It turns out that the data 
cannot be explained in NLO by using px as scale, unlike the case of inclusive charged- 
hadron production. Both choices of scales fi = M = Mf = p^ and 2pT give cross sections 
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Figure 12: Scale dependence of the charged-hadron production cross section for ^/s = 630 GeV, 
\y\ < 1.0 and pT = 5,10 and 30 GeV. The sohd hues show the scale dependence of the NLO 
cross section. The dotted and dashed lines show the behaviour of the bona fide LO cross section 
and the LO one with NLO structure functions and as, respectively. 
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Figure 13: Rapidity distribution of the charged-hadron production cross section at ^/s = 
630 GeV, scales fixed to pr, for px = 5,10 and 30 GeV. The upper curves show the y shape 
of the NLO cross section, the lower curves the shape of the LO one with structure functions and 
as at the NLO level. 
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Figure 14: Inclusive tt^ production cross section for ^/s = 540 GeV averaged over the rapidity 
range 1.0 < y < 1.8 compared to the experimental data obtained by the UA2 Collaboration. 
The short-dashed, solid and long-dashed lines correspond to the full NLO predictions for scales 
H, M and Mf all set equal to pT, 2pT and ApT, respectively. 
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definitely above the data. A reasonably good fit can be obtained only if we increase the 
value of these scales up to 4pT- On the other hand, the fragmentation functions into 
7r°'s might need to be revised, in particular the gluon fragmentation function, which is 
not much constrained by the e~^e~ data used in More work on the 7r° fragmentation 
functions is needed before a final conclusion concerning the scale appropriate to the vr*^ 
production cross section can be drawn. 
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5 Summary and Conclusions 



We have calculated inclusive single-charged-hadron and single-vr" cross sections in NLO 
and compared our results with experimental data from the UA2 and CDF Collaborations. 
The overall agreement is satisfactory, surprisingly, also in the small-px region. At large p^, 
where QCD could be tested cleanly, experimental data are unfortunately rare and plagued 
by large errors. More accurate measurements in this region would obviously provide a 
better probe of the QCD predictions. 

As expected, the scale dependence in NLO is reduced with respect to the one present 
in LO. Therefore, within the scale variation in NLO, our results can be considered as an 
absolute prediction for single-charged- and single-neutral-hadron production in pp colli- 
sions. The scale which yields the best agreement of the inclusive charged-hadron produc- 
tion cross section to the experimental data is equal to px- Larger values, up to 4pT, are 
needed to fit the data for inclusive 7r° production. 

A drawback of our approach is that we still use fragmentation functions in LO. While 
this will be improved in the future when NLO fragmentation functions become available, 
we observe that already now the scale dependence of our results is reduced with respect 
to the one present in the LO calculation. We deduce that some compensation concerning 
the factorization scale Mf already works to a large extent. 

When comparing our results to the LO ones with LO structure functions and one-loop 
as, we find that they roughly agree at scales just below px- We consider this agreement 
as accidental, since they deviate from each other at different scales. The fact that the 
NLO hard scattering contributions are very important can be clearly seen when the NLO 
results are compared with the LO ones with as and structure functions taken in NLO. In 
this case, for scales equal to p^, we find a K factor approximately 1.5. 

We plan to extend this NLO approach to inclusive hadron production in jp and ep 
collisions. We expect in this case a somewhat different pattern due to the presence of the 
direct component in the photon structure function. Further work is also needed in the 
improvement of the fragmentation functions where a thorough analysis, similar to the one 
performed for the structure functions, is called for. 
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